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We have developed a multi-objective optimization (MOO) procedure to construct modified- 
embedded-atom-method (MEAM) potentials with minimal manual fitting. This procedure has been 
applied successfully to develop a new MEAM potential for magnesium. The MOO procedure is 
designed to optimally reproduce multiple target values that consist of important materials proper- 
ties obtained from experiments and first-principles calculations based on density-functional theory 
(DFT). The optimized target quantities include elastic constants, cohesive energies, surface ener- 
gies, vacancy formation energies, and the forces on atoms in a variety of structures. The accuracy of 
the new potential is assessed by computing several material properties of Mg and comparing them 
with those obtained from other potentials previously published. We found that the present MEAM 
potential yields a significantly better overall agreement with DFT calculations and experiments. 

PACS numbers: 34.20.Cf, 61.43.Bn, 61.72.Ji, 62.20.Dc, 68.35.-p 



I. INTRODUCTION 

Molecular dynamics simulations are effective tools used 
to study many interesting phenomena such as the melting 
and coalescence of nanoparticles at the atomic scaleJ^ 
These atomistic simulations require accurate interaction 
potentials to compute the total energy of the system, and 
first-principles calculations can provide the most reliable 
interatomic potentials. However, realistic molecular dy- 
namics simulations often require an impractical number 
of atoms that either demands too much computer mem- 
ory or takes too long to be completed in a reasonable 
amount of time. One alternative is to use empirical or 
semi-empirical interaction potentials that can be evalu- 
ated efficiently. 

The modified-embedded-atom method (MEAM) pro- 
posed by Baskes et alj'S.iiJ, is one of the most widely 
used methods using semi-empirical atomic potentials to 
date. The MEAM is an extension of the embedded- 
atom method (EAM) to include angular forcesi^!^ The 
MEAM and EAM use a single formalism to generate 
semi-empirical potentials that have been successfully ap- 
plied to a large variety of materials including fee, bcc, 
hep, diamond-structured materials and even gaseous el- 
ements, to produce simulations in good agreement with 
experiments or first-principles calculations ^^^i^ 

Despite its remarkable successes, one of the most no- 
table difficulties in using MEAM is that the construction 
of the MEAM potentials involves a lot of manual and ad 
hoc fittings. Because of the complex relationship between 
the sixteen MEAM parameters and the resultant behav- 
ior of a MEAM potential, a traditional procedure for con- 
structing a MEAM potential involves a two-step iterative 
process. First, a single crystal structure, designated as 



the reference structure, is chosen and the MEAM pa- 
rameters are fitted to construct a MEAM potential that 
reproduces a handful of critical materials properties of 
the element in the reference structure. Second, the new 
potential is tested for its accuracy and transferrability by 
applying it to atoms under circumstances not used during 
its construction phase. These systems include different 
crystal structures, surfaces, stacking faults, and point de- 
fects. If the validation is not satisfactory, one needs to 
go back to the first step and adjust the parameters in 
a way that improves the overall quality of the potential. 
Although this iterative method does work eventually in 
many cases, it is a very tedious and time-consuming. Er- 
colessi and Adams overcame this shortcoming for EAM 
potentials by developing a force-matching method that 
fits the EAM potential to ah initio atomic forces of many 
atomic configurations including surfaces, clusters, liq- 
uids and crystals at finite temperaturei^ Later, the force- 
matching method was extended to include many other 
materials properties such as cohesive energy, lattice con- 
stants, stacking fault energies, and elastic constants i^^d^ 
Furthermore, several different MEAM potentials for the 
same element often develop and an objective and quan- 
titative method to measure the relative quality of each 
potential would be helpful for the researchers who want 
to choose one of these potentials. 

In this work, we extend the force-matching method 
to develop a multi-objective optimization (MOO) proce- 
dure to construct MEAM potentials. Most realistic op- 
timization problems, particularly in engineering, require 
the simultaneous optimization of more than one objec- 
tive function. For example, aircraft design requires si- 
multaneous optimization of fuel efficiency, payload and 
weight calls for a MOO procedure. In most cases, it 
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is unlikely that the different objectives would be opti- 
mized by the same parameter choices. Therefore, some 
trade-off between the objectives is needed to ensure a 
satisfactory design. Stadler-'^^ introduced the concept of 
Pareto optimality^'^ to the fields of engineering and sci- 
ence. The most widely used method for multi-objective 
optimization is the weighted sum method. A comprehen- 
sive overview and comparison of different MOO methods 
can be found in Ref . 

The composite objective function also provides an un- 
biased measure to quantify the relative quality of differ- 
ent MEAM potentials. We apply the procedure to de- 
velop a new MEAM potential for magnesium. The new 
Mg MEAM potential will be compared with previously 
published Mg potentials. 

We chose Mg bacause of its increased importanta- 
nce in many technological areas, including the aerospace 
and automotive industries. Due to the lower mass 
densities of magnesium alloys compared with steel and 
aluminum and higher temperature capabilities and im- 
proved crash-worthiness than plastics, the use of magne- 
sium die castings is increasing rapidly in the automotive 
industry^ii^iii 

Empirical potentials for Mg have been previously pro- 
posed by several groups. In 1988, Oh and Johnson devel- 
oped analytical EAM potentials for hep metals such as 
MgfiS, Igarashi, Kanta and Vitek— (IKV) also developed 
interatomic potentials for eight hep metals including Mg 
using the Finnis-Sinclair type many-body potentials.— 
Pasianot and Savinc^^ proposed improved EAM poten- 
tials for Mg based on IKV's fitting scheme. Baskes and 
Johnson* have extended the modified embedded atom 
method (MEAM)^>^i^ to hep crystal structures. Later, 
Jelinek et al. improved this potential as a part of the 
MEAM potentials for Mg-Al alloy system.— Liu et al. 
used the force-matching method to develop an EAM po- 
tential for Mg.iS 

The paper is organized in the following manner. In 
Sec. im we give a brief review of the MEAM. In Sec. [TTIl 
the procedure for determination of the MEAM param- 
eters is presented in detail. In Sec. IIVI we assess the 
accuracy and transferability of our MEAM potential and 
make comparisons to other previously published poten- 
tials. 



The energy of atom i consists of the embedding energy 
and the pair potential terms: 

E,^F,ip,) + ^J2c^,,{n,). (2) 

Fi is the embedding function of atom i; pi is the back- 
ground electron density at the site of atom i; and (rij) 
is the pair potential between atoms i and j separated by 
a distance rij. The embedding energy Fi(pi) represents 
the energy cost to insert atom i at a site where the back- 
ground electron density is pi. The embedding energy is 
given in the form 

F,{p,)^AE^P,Hp,), (3) 

where the parameters and Ai depend on the element 
type of atom i. The background electron density pi is 
given by 

(0) 

Pi = ^G(r,), (4) 

where 

k=l \pi ) 

and 

G(r) = ^/TTT. (6) 

The zeroth and higher order densities, p^i \ pf^ i ^^'^ 
pf^ are given in Eq. ([5]). The composition-dependent 
electron density scaling p° is given by 

p^ = p,,Z,,G{Tf), (7) 



II. METHODOLOGY 



MEAM 



The total energy i? of a system of atoms in the 
MEAM2^ is approximated as the sum of the atomic en- 



ergies 



E = Y,E, 



(1) 



where p^o is an element-dependent density scaling, Z^o 
is the first nearest-neighbor coordination of the reference 
system, and F ■'^^ is given by 



prcf _ J_\^ji.k) (k) 
"^^ fc=l 



(8) 



where s"^^ is the shape factor that depends on the ref- 
erence structure for atom i. Shape factors for various 
structures are specified in the work of Baskes*. The par- 
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tial electron densities are given by 



(9a) 
(9b) 
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(9c) 



Hp. 



a(3) '^ijoi'^ij 



-Is 



E a(3)^;Xja o 



(9d) 



where Ty o, is the a component of the displacement vector 
from atom i to atom j. Sij is the screening function 
between atoms i and j and is defined in Eqs. (fT6|) . The 
atomic electron densities are computed as 



Pi^ V'l]) = P*oexp 



(10) 



where t\ is the nearest-neighbor distance in the single- 
element reference structure and (if^ are element- 
dependent parameters. Finally, the average weighting 
factors are given by 



r(fe) _ 1 ,(fc) a(0)„ 

~ (0) '^i 



(11) 



where tj*^^ is an element-dependent parameter. 
The pair potential is given by 

1 



2EUn,)-F., 




(12) 

(13) 

(14) 
(15) 



where a^- is an element-dependent parameter. The subli- 
mation energy Efj , the equilibrium nearest-neighbor dis- 
tance r^j , and the number of nearest-neighbors Z^j are 
obtained from the reference structure. 



The screening function Sij is designed so that Sij = 1 
if atoms i and j are unscreened and within the cutoff 
radius Vc, Sij = if they are completely screened or 
outside the cutoff radius, and varies smoothly between 
and 1 for partial screening. The total screening function 
is the product of a radial cutoff function and three-body 
terms involving all other atoms in the system: 



(16a) 
(16b) 

(16c) 
(16d) 

(16e) 
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Sikj — /( 



ikj 



1 + 2 
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2 2 

'^ik ~ '^jk 



fcix) 



1 X>1 

[1- (l-a;)4)]' 0<a;<l 
x<0 



Note that Cmin and Cmax can be defined separately for 
each i-j-k triplet, based on their element types. The 
parameter Ar controls the distance over which the radial 
cutoff function changes from 1 to near r = rc- 



B. Multi-objective Optimization 

A generic multi-objective optimization (MOO) prob- 
lem can be formulated as2i2^: 



min J(x) s.t. x € S 
where J — [Jii^) ■ ■ ■ Jm(x)]'^ 
X = [xi • • ■x„]'^ 



(17) 



Here, J is a column vector of m objectives, whereby 
Ji (z m. The individual objectives are dependent on a 
vector X of n design variables in the feasible domain S. 
The design variables are assumed to be continuous and 
vary independently. Typically, the feasible design domain 
is defined by the design constraints and the bounds on the 
design variables. The problem is to minimize all elements 
of the objective vector simultaneously. The most widely 
used method for MOO is scalarization using the weighted 
sum method. The method transforms the multiple objec- 
tives into an aggregated scalar objective function J that 
is the sum of each objective function Ji multiplied by a 
positive weighting factor Wi: 



J(x) 



W,;Ji(x). 



(18) 



In this work, the overall goal is to develop a MEAM 
potential for Mg. The individual objective functions are 
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constructed from the normalized differences between the 
MEAM-generated values and the target values: 



J.(x) 



(19) 



Here, Qi is the physical quantity computed using the 
current MEAM potential parameters and is the tar- 
get value to reproduce. The target values are usually 
experimental values, but the computed values from the 
first-principles method are chosen when the experimen- 
tal data are not available. The normalization factor Q* 
is a typical value for the given materials parameter and 
often Q* = Qi- The overall objective function J(x) can 
be minimized using usual multi-dimensional otimization 
routines. To avoid unnecessary complications, we used 
the downhill simplex method,^^ which requires only func- 
tion evaluations, not derivatives. 



a crystal structure: 



^'coh — 



N 



(20) 



where Stot is the total energy of the system, N is the 
number of Mg atoms in the system, and i?atom is the 
total energy of an isolated Mg atom. The cohesive ener- 
gies of Mg atoms in hep, fee, and bcc crystal structures 
for several atomic volumes near the equilibrium atomic 
volumes were calculated. Fig. [1] is an example of the 
cohesive energy plot of Mg atoms as a function of the 
lattice constant. The minimum of this curve determines 
the equilibrium lattice constant gq and equilibrium cohe- 
sive energies E'hcp = ^-hcp in Table HIl 



2. Elastic constants 



III. POTENTIAL CONSTRUCTION 
PROCEDURE 

We used the MOO procedure to develop a new set of 
MEAM parameters that improves the overall agreement 
of MEAM results with experiments or ab initio calcula- 
tions. Our previously published MEAM parameters for 
Mg2^ served as the basis for the present work. 

All ab initio total-energy calculations and geome- 
try optimizations are performed within density func- 
tional theory (DFT) using ultrasoft pseudopotentials 
(USPP)^'' as implemented by Kresse et. a L^^'^^ For 
the treatment of electron exchange and correlation, we 
use local-density approximation (LDA)^i2i. The Kohn- 
Sham equations are solved using a preconditioned band- 
by-band conjugate-gradient (CG) minimization,— The 
plane-wave cutoff energy is set to at least 300 cV in all 
calculations. Geometry relaxations are performed until 
the energy difference between two successive ionic opti- 
mizations is less than 0.001 eV. The Brillouin zone is sam- 
pled using the Monkhorst-Pack scheme^ and a Fermi- 
level smearing of 0.2 eV was applied using the Methfessel- 
Paxton methodj^i 

The objectives used in this work include equilibrium 
hep lattice constants a and c at K, the cohesive energy, 
elastic constants, vacancy formation energy, surface en- 
ergies, stacking fault energies, and adsorption energies. 
We also used the forces on Mg atoms in structures equi- 
libriated at six different temperatures. The final MEAM 
parameters obtained from the MOO procedure are listed 
in Table m Table HIl shows the complete list of objectives 
optimized to construct the MEAM potential parameters 
for Mg and their weights. 

1. Cohesive energies 

The cohesive energy of Mg atom is defined as the heat 
of formation per atom when Mg atoms are assembled into 



Hexagonal crystals have five independent elastic con- 
stants: Cii, Ci2, Ci3, C33, and €44^ The elastic 
constants are calculated numerically by applying small 
strains to the lattice. For small deformations, the rela- 
tionship between deformation strain and elastic energy 
increase in an hep crystal is quadraticf^ 

1. AU = (5^(Cii + C12), for deformation x' = a; + (5 • a; , 
1/ = y + d -y, 

2. AU = (5^(Cii — C12), for deformation x' = a; -f (5 • x, 

y' = y-s ■ y, 

3. AU = 3^0^3/2, for deformation z' = z + S ■ z, 

4. AU = (52(2Cii + C33 + 2Ci2 + 4Ci3)/2, for defor- 
mation x' = X + S ■ X, y' = y + 6 ■ y, z' = z + S ■ z, 

5. AU — (5^6*44/2, for deformation z' = z + S ■ x, 

where unprimcd (primed) are the coordinates of the lat- 
tice before (after) deformation. AU is the elastic energy 
due to the deformation, and S is the small strain applied 
to the lattice. We follow the procedure described by Mehl 
et al."^^ and apply several different strains ranging from 
—2.0% to -1-2.0%. The elastic constants are obtained by 
fitting the resultant curves to quadratic functions. We 
found that this method gives much more stable results 
than using one strain valuei^,. 

3. Surface formation energies 

Surface formation energy per unit surface area is de- 
fined as 

7 = [Etot - Ne) /A, (21) 

where i?tot is the total energy of the system with a sur- 
face, N is the number of atoms in the system, e is the 
total energy per atom in the bulk, and A is the sur- 
face area. Table HI] lists the surface formation energies 
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used in this study. The (1010) surface of hep crystals 
can be terminated in two ways, either with a short first 
interlayer distance di2 ("short termination") or with a 
long di2 ("long termination") (See, for example. Fig. 2 
of Ref. l37r ). In this study, we only included the results 
for the short terminated surface, since it is known to 
be enegetically more favorable over the long terminated 
surfacc'^^ in agreement with experimental observations in 
Be(lOTO) and other hep metalsi^ 

4- Stacking fault energies 



crystal structure. The positions of atoms are randomly 
disturbed from their equilibrium positions and 10 000 
steps of molecular-dynamics (MD) simulations with a 
timestep of At — 2.5 ps were performed to equilibrate 
each structure for different temperatures. In each MD 
run, we used Mg MEAM potential by Jelinek et alM If 
no MEAM potential were available for MD simulations, 
one could use an intermediate MEAM potential that is 
generated with this MOO procedure without the force 
data. The potential should be adequate enough to ob- 
tain a reasonable set of structures. 



by 



Stacking fault formation energy per unit area is defined 



Esf = {Etot - Ne) /A, 



(22) 



where Etot is the total energy of the structure with a 
stacking fault, N is the number of atoms in the sys- 
tem, £ is the total energy per atom in the bulk, and 
A is the unit cell area that is perpendicular to the 
stacking fault. For Mg, four stacking fault types from 
the calculation of Chetty et al.^^ were examined. The 
sequences of the atomic layers within the unit cell of 
our simulations are: h = ABABABCBCBCB, h = 
ABABABCACACB, T2 = ABABABCBABAB, and 
E = ABABABCABABAB. We note that the unit cells 
for Ii and I2 contain two stacking faults and the quan- 
tities obtained from Eq. ((22| must be divided by two to 
obtain the correct formation energies. 



5. Vacancy formation energies 

The formation energy of a single vacancy E^, 
fined as the energy cost to create a vacancy: 



E. 



Etot[N]-Ne, 



is de- 



(23) 



where £^tot[-^] is the total energy of a system with N 
atoms containing a vacancy, and e is the energy per atom 
in the bulk. 



Atomic Forces 



For forces, the objective functions are defined as: 



J.(x) 



(((F-F")2) 
(((F°)2))V2 



,1/2 



(24) 



where F are the force vectors on atoms calculated using 
the MEAM while F° are the force vectors from DFT 
method. (((F°)^))^/^ represents the root-mean-square of 
the DFT force, and {{{¥ - Y^f))^/"^ is the root-mean- 
square of the error in the force. 

To obtain the force data, initial atomic structures that 
contain 180 Mg atoms were created from the bulk hep 



IV. RESULTS AND DISCUSSION 

The hep structure was chosen as the reference structure 
for Mg. The final MEAM parameters obtained from the 
MOO procedure are listed in Table |T1 



A. Materials properties 

Table |lT] lists various materials properties of Mg se- 
lected as the objectives to be optimized in constructing 
Mg MEAM potential, along with experimental data and 
ah initio data. It also shows how well each objective has 
been optimized. Results from other previously published 
Mg potentials are also listed in the table for comparison. 
Table mi also shows the weight of individual objectives Wi 
chosen to optimize the present potential. The underlined 
quantities are the target values chosen for the MOO pro- 
cedure. Whenever possible, the experimental values are 
chosen as the target values. If the experimental values 
are not available or known to be unreliable, the computed 
values from the first-principles method are used. 

The present MEAM potential reproduces the experi- 
mental lattice constant, the c/a ratio, and the cohesive 
energy near perfectly. Fig. [T] shows the cohesive en- 
ergies of Mg atoms in hep crystal structure compared 
with those obtained from the Rose universal equation of 
state'*^ based on the experimental lattice constant, cohe- 
sive energy and bulk modulus. It shows a good agreement 
between the two sets of data. We also note that the se- 
quence of the structures is predicted correctly in the order 
of stability by the present Mg MEAM potential as shown 
in Table ini 

The surface formation energies of the two common 
low-index surfaces of hep Mg crystals are in good agree- 
ment with the experimental values, representing a signif- 
icant improvement over the previously published MEAM 
potentials 

As pointed out by Liu et alJ^, the stacking fault ener- 
gies are difficult quantities for an emprirical potential to 
reproduce because they only depend on long range inter- 
actions beyond second nearest-neighbor distances in hep 
crystals. The present MEAM potential shows a substan- 
tial improvement over the previously published MEAM 
potential by Hu et al4^ The stacking fault energies are 
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TABLE I: Set of the MEAM potential parameters for Mg. See Sec. Ill Al for the definition of these parameters and their usage. 
The hep structure was chosen as the reference structure for Mg. 

£°[eV] ao[A] A a /j'^) /j'^) /j'^) ^(o) ^(i) ^(2) ^(3) ^^^^ ^^.^ 

1.51 3.20 1.14 5.69 2.66 -0.003 0.348 3.32 1.00 8.07 4.16 -2.02 3.22 1.10 5.0 0.353 



TABLE IL The objectives optimized to construct the MEAM potential parameters for Mg along with the target values and 
their weights. Comparisons to other Mg potentials are also made. Ecoh is the cohesive energy, B is the bulk modulus, 7 is the 
surface energy, E^i is the stacking fault formation energy, -Evac is the relaxed vacancy formation energy, and AF is the relative 
error in forces. The underlined quantities are the target values chosen for the MOO procedure. 
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consistently underestimated by the present MEAM po- 
tential compared to the results of the DFT calculations, 
while the results by the EAM potential from Ref. [l^ 
are consistently overestimated. Table |TT] also shows that 
the formation energy of single vacancies from DFT cal- 
culation is reproduced quite reasonably by the present 
MEAM potential. 

Table |TT] also shows the force-matching against the ab 
initio forces database. It shows that the MEAM po- 
tential from the present work reproduces more accu- 
rate forces on atoms compared to the previous MEAM 
potential—. 

B. Additional materials properties 

To validate the present MEAM potential further, we 
calculated a few additional materials properties of Mg 



TABLE 111: The additional materials properties of Mg that 
are not used as objectives for the constructuction of the po- 
tential. Comparisons are made to other Mg potentials and 
experiments. -Bads is the adsorption energy; BJ"' is the for- 
mation energies of interstitial point defects. All energy values 
are given in eV. 



Property 


DFT 


MEAM" 


Jelinek' 


Bads [(0001)] 


-0.81 


-1.46 


-1.50 


Sads [(1010)1 


-1.21 


-1.52 


-1.56 


Bj"' (octahedral) 


2.36 


1.20 


1.29 


Bj"' (tetrahedral) 


2.35 


1.41 


1.53 



"MEAM potential from the present work 
''MEAM potential from Ref. [H 



that were not used as objectives during the construction 
of the potential. We obtained the adsorption energies of 
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Lattice Constant (A) 



FIG. 1: The cohesive energies as a function of the lattice 
constant a for Mg atoms in hep crystal structure compared 
with the ones obtained from the Rose equation. The data 
points are computed with the present MEAM potential while 
the curve is obtained from the Rose equation. 

a single Mg atom on different surfaces and the formation 
energies of interstitial defects as listed in Table Hill 
The adsorption energy of a single adatom i?ads is given 

by 

^-ads — Etot ~ ^-surf — -E^atom , (25) 

where Etot is the total energy of the structure with the 
adatom adsorbed on the surface, i?surf is the total en- 
ergy of the surface without the adatom, and i?atom is the 
total energy of an isolated atom. On both (0001) and 
(1010) surfaces, we placed a single Mg atom at the site 
where the atoms of the next layer would normally sit. 
The structures were then relaxed to determine the ad- 
sorption energies. Table IIIII shows that the adsorption 
energies on two Mg surfaces are quite well reproduced 
by the present MEAM potential. The present Mg po- 
tential gives slightly better adsorption energies than the 
previously published MEAM potential^^. 

The formation energy of an interstitial point defect 
E"™* is given by 

E't^ = ^tot [A^ + 1] - {Etot m + e), (26) 

where E'tot [A] is the total energy of a system with A^ Mg 
atoms, i<^tot[A -I- 1] is the total energy of a system with 
A' atoms plus one Mg atom inserted at one of the inter- 
stitial sites, and e is the total energy per Mg atom in its 
most stable bulk structure. Interstitial atom formation 
energies were calculated for Mg at octahedral and tetra- 
hedral sites. Atomic position and volume relaxation were 
performed. The results of these calculations are listed in 
Table Hill to be compared with the results from the DFT 
calculations. The present MEAM potential predicts cor- 
rect signs for these energies although the magnitudes are 
about half of those predicted by DFT. MEAM potentials 
predict that the octahedral site will be more stable than 
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FIG. 2: The internal energies of Mg crystal in hep structure 
as a function of temperature. The energies are obtained from 
the ensemble average of the MD simulations of five structures 
containing 448 Mg atoms. 

the tetrahedral site, while the DFT calculations indicate 
that both sites will have nearly the same formation en- 
ergies. 

C. Thermal properties 

To validate the new potential for molecular dynamics 
simulations, we calculated the melting temperatures of 
pure Mg crystals. We followed a single-phase method 
as described by Kim and Tomanek," in which the tem- 
perature is increased at a constant rate and the internal 
energy of the system is monitored. Fig. [2| shows the in- 
ternal energies of Mg crystal in hep structure as a func- 
tion of temperature. The plot was obtained from the en- 
semble average of five hep structures containing 448 Mg 
atoms. The initial velocity vectors were set randomly 
according to the Maxwell-Boltzmann velocity distribu- 
tion at r = 100 K. The temperature of the system was 
controlled by using a Nose-Hoover thermostat 4^1^ It is 
clearly seen from Fig. [2| that the internal energy curve 
makes an abrupt transition from one linear region to an- 
other, marking the melting point. Using this method, 
we obtained 920 K as the melting temperature of Mg 
crystals. This result is in good agreement with the ex- 
perimental value of 923 K. Our result represents a sub- 
stantial improvement in accuracy from 745 K obtained 
from a previously published EAM potential^° or 780 K 
from a MEAM potential.^ 

V. CONCLUSIONS 

In this study, we developed a multi-objective optimiza- 
tion procedure to construct MEAM potentials with min- 
imal manual fitting. We successfully applied this proce- 
dure to develop a set of MEAM parameters for Mg in- 
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teratomic potential based on first-principles calculations 
within DFT. The validity and transferability of the new 
MEAM potentials were tested rigorously by calculating 
the physical properties of the Mg systems in many dif- 
ferent atomic arrangements such as bulk, surface, point 
defect structures, and molecular dynamics simulations. 
The new MEAM potential shows a significant improve- 
ment over previously published potentials, especially for 
the atomic forces and melting temperature calculations. 
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